LHP1-mediated epigenetic buffering of subgenome diversity and defense responses confers genome plasticity and adaptability in allopolyploid wheat

Polyploidization is a major driver of genome diversification and environmental adaptation. However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity. By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenome-diversified genes. Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis. Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression. The simultaneous inactivation of LHP1 homoeologs by CRISPR–Cas9 confers robust stripe rust resistance in wheat seedlings. The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution. These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat. Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.

Polyploidization is a major driver of genome diversification and environmental adaptation.However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity.By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenomediversified genes.Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis.Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression.The simultaneous inactivation of LHP1 homoeologs by CRISPR-Cas9 confers robust stripe rust resistance in wheat seedlings.The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution.These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat.Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.
Polyploidization, which is a recurring event during evolution, is a major driver of genome diversification that promotes adaptive evolution [1][2][3] .It is pervasive in both plants and animals.Many eukaryotic genomes have a polyploid ancestry, and a substantial proportion of crops underwent recent polyploidization events 4,5 .The rise of allohexaploid wheat (Triticum aestivum, 2n = 6x = AABBDD) approximately 10,000 years ago shaped the evolution of modern humans 6 .However, the convergence of different genomes may not necessarily result in heterosis in polyploids.Instead, the genetic heterogeneity between subgenomes (hereafter referred to as 'subgenome diversity') may result in genetic conflicts.Typical examples include competition between parental genomes and the rapid loss or repression of homoeologous gene copies (biased fractionation or repression) [7][8][9][10] .Persistent aneuploidy and decreased fertility are generally associated with nascent allohexaploid plants 11,12 .How sequence diversity in polyploids is effectively interpreted and regulated to enhance environmental plasticity is an important but under-explored issue.
Genetic diversity and environmental stimuli respectively represent internal and external stresses affecting organisms.Phenotypes need to be robust in response to internal and external changes, which requires a buffering system to ensure developmental flexibility ('canalization') 13,14 .Exposure to long-term environmental stress can reveal genetic variants, which may be achieved by repressing buffering systems.These uncovered variants mediating adaptive phenotypes are selected and genetically fixed (assimilation) 15 .In eukaryotes, chromatin and the associated epigenetic mechanisms represent a typical buffering system that stabilizes transcription and cellular homeostasis against internal and external changes, while also reprogramming the transcriptome in response to developmental or environmental cues.Incorporating diverse genomes into a single nucleus yields large-scale epigenomic changes.Epigenetic diversity across subgenomes and the flexible interaction with the transcriptional machinery may provide polyploid wheat with a selective advantage [16][17][18][19][20] .Recent advances detected wide-spread subgenome-unbalanced transcription 21,22 , which is mediated by the interplay between genetic and epigenetic diversity, including open chromatin 23,24 , histone modifications 25,26 , and specific transcription factor (TF) binding densities 27 , which may contribute to phenotypic adaptation.However, the available evidence is mostly based on statistical associations between high-throughput data, and a causal relationship between epigenetic programming and phenotypic plasticity remains to be determined.
In addition to influencing transcriptional diversity 17,18 , recent investigations on Arabidopsis 28 and fungal populations 29,30 indicated that epigenetic architecture directly affects sequence variation and adaptation rates.This effect might be conserved in other eukaryotes, potentially influencing the pace and nature of evolution.Identifying the factors mediating genetic and epigenetic interactions is critical for characterizing the mechanisms by which subgenome diversity and epigenetic plasticity cooperatively promote environmental adaptation and for clarifying and harnessing polyploid plasticity.
In this study, by analyzing subgenomic variations and regulomics data as well as performing genetic screening and functional validation experiments, we revealed how the epigenetic system mediates genotype-environment interactions that determine polyploid wheat phenotypes.The buffering effect of this system on subgenome diversity contributes to developmental plasticity in response to external stimuli and promotes the mutation-selection-fixation cycle of adaptive evolution.

LHP1 is a master repressor of subgenome-diversified genes
Polyploid evolution depends on subgenome diversity 31,32 .Identifying the factors regulating subgenome variations is vital for elucidating the causal mechanism underlying polyploid adaptation and evolution.Among the 107,891 high-confidence full-length genome-wide gene models in common wheat, approximately one-third were single copy (sc)-triads (1:1:1 correspondence across subgenomes) (Fig. 1a), whereas the rest were non-sc-homoeologs (varying homoeolog numbers).To reveal the main regulators of non-sc-homoeologs, we collected and generated chromatin immunoprecipitation sequencing (ChIP-seq) and DNA affinity purification sequencing (DAP-seq) data for 191 genomewide DNA-or chromosome-binding factors and then searched for factors that preferentially regulate non-sc-homoeologs (Supplementary Data 1).On the basis of LASSO regression 33 , which is a supervised regularization method used in machine learning because of its powerful built-in feature-selection capability (Fig. 1b), LHP1 was identified as the primary trans-factor regulating non-sc-homoeologs (Fig. 1c, top), especially the homoeologs absent in one subgenome and multicopies in other subgenome(s), including N:0:N, 0:N:N, 0:1:N (N indicates a minimum of one additional paralog per respective subgenome) (Fig. 1d).Increases in the copy number coincided with increases in the degree of LHP1 binding (Supplementary Fig. 1).The AP2-type TFs were predicted to be among the major regulators of sc-triads (Fig. 1c, bottom).Consistent with these analyses, a large proportion of these enriched homoeolog groups were targeted by LHP1, whereas only 6.8% sc-triads are occupied by LHP1 (Fig. 1e).
In plants, LHP1 is a member of the Polycomb group proteins (PcGs) 34 and is responsible for recognizing and spreading the repressive H3K27me3 mark [35][36][37][38] .The specific recruitment of LHP1 is critical for plant development [39][40][41][42][43] .We fine-mapped LHP1-binding sites throughout the common wheat genome on the basis of a ChIP-seq analysis.These binding sites were typically co-localized with the H3K27me3 mark (Fig. 1f), as illustrated by the genomic tracks around VERNALIZATION1 (VRN1) (Fig. 1g).A phylostratigraphic analysis 44 , which clustered the genes of common wheat based on their orthologous relationships with nine other species (see "Methods"), revealed an apparent burst of H3K27me3 sites in Triticeae (Fig. 1h), suggestive of an increase in LHP1-mediated epigenetic regulation in Triticeae.
To functionally validate LHP1 and explore the mechanism by which it controls subgenome-diversified genes, we inactivated all three LHP1 homoeologs in the common wheat variety 'JW1' using the CRISPR-Cas9 editing system and performed ChIP-seq and RNA-seq analyses to assess the effect of a loss-of-function mutation to LHP1 on the expression and epigenetic status of subgenome-diversified genes.Two homozygous triple mutant (Talhp1-abd) lines with different frame-shift mutations in all three homoeologs were developed (Fig. 2a and Supplementary Fig. 2).There were apparently fewer H3K27me3 marks in the triple mutant than in the single and double mutants (Supplementary Fig. 3).Similar to the corresponding Arabidopsis mutants, the mutants lacking LHP1 in all three subgenomes flowered early and exhibited dwarfism at maturity (Fig. 2b).We quantified the H3K27me3 levels in wild-type and mutant plants using a modified spike-in method that involved the addition of the same amount of Arabidopsis chromatin before ChIP to normalize samples 45 (details are provided in the "Methods" section).The analysis of the Talhp1-abd mutant detected a global decrease in H3K27me3 surrounding genes, most of which were non-sc-homoeologs (Fig. 2c, d).The comparison with the transcriptomic changes in Talhp1abd revealed a close association between the decreased H3K27me3 levels and the increased target gene expression levels (Fig. 2e).Considered together, these observations showed that LHP1, through its capacity to control the H3K27me3 level, is a master repressor of subgenome-diversified gene activity.

LHP1-mediated H3K27me3 is involved in maintaining and promoting genetic diversity across subgenomes
There is accumulating evidence that the evolutionary rate of a gene is predominantly influenced by its expression level rather than its functional importance 46 .The negative correlation between gene expression and the evolutionary rate exists in all three domains of life 46 .The strength and breadth of the effects of LHP1 on the repression of subgenome-diversified gene activity suggest that LHP1 may influence evolutionary processes.We examined the relationship between LHP1 binding and target sequence variations at the population level.
The subgenome diversity in common wheat has two major sources: captured from diploid progenitors (i.e., variations already present across diploids) and newly generated in polyploid populations 31,32 .We first compared the extent of the captured variations (reflected by the differences between subgenomes of the same individual) and polyploidy-generated diversity (reflected by the population diversity π of the same subgenome; i.e., genetic diversity across hexaploid individuals 47 ) in a pairwise framework (Fig. 3a).There is a high correlation between captured diversity and polyploidy-generated diversity.The variable genome regions in progenitors tended to further diverge after allopolyploidization (region 2 in Fig. 3a), representing selective neutral or weakly deleterious loci likely due to genetic drift.
Next, the enrichment of LHP1 target loci in the above pairwise comparisons was determined.The LHP1-binding sites were highly enriched in regions with high levels of captured diversity, which were further diversified in polyploid populations (Fig. 3b).Among the various epigenetic marks, H3K27me3 was significantly over-represented in these loci (Fig. 3c).Accordingly, LHP1-mediated H3K27me3 repression appears to be closely associated with an increase in subgenomic diversity at the population level, possibly because it promotes neutral   variations and protects the diversified genes from strong purifying selection 46 .The genomic tracks in Fig. 3d reflect the quantitative associations among LHP1 binding, H3K27me3 modifications, and population diversity.In accordance with these findings, populationlevel investigations of defective epigenetic enzymes in fungi demonstrated epigenetic repression, typically via H3K27me3, promotes mutations and adaptations 29,30 .Moreover, a recent genomewide characterization of crossover-active regions and relevant factors in common wheat suggested that H3K27me3 is positively associated with crossovers 48 .These findings and evidence imply that LHP1  :1:1) and non-sc-homoeologs with varying number of homologs across subgenomes in common wheat.b Workflow for predicting the core factors regulating sc-triads (1:1:1) and non-sc-homoeologs across subgenomes.c LASSO coefficients of the trans-factor binding-related regulation of non-sc-homoeologs (top) and sc-triad (bottom) genes.The coefficient value represents the relative importance of the trans-factor contributing to regulation.d Enrichment of different homoeologous groups among LHP1-targeted genes.e Fractions of homoeolog groups targeted by LHP1.f ChIP-seq read densities of LHP1 and H3K27me3 surrounding H3K27me3 peaks.g Genomic tracks illustrating the co-occupation of LHP1 and H3K27me3 ChIP-seq signals surrounding VRN1-5A and its homoeologs.h Enrichment of the typical epigenetic marks in eight phylostrata corresponding to the phylogenetic internodes, which are bordered by vertical grids and denote sets of genes whose founder genes originated in the corresponding evolutionary periods (e.g., the left-most phylostratum represents the common ancestor of dicots and monocots).The color represents the enrichment score.High enrichment scores reflect the expansion of the targets with a given epigenetic mark in the corresponding phylostrata.A value = (JW1+lhp1-abd)/2 590 genes 13,330 genes repression protects and potentially promotes the accumulation of neutral or weakly deleterious variations across subgenomes.

Direct repression of subgenome-diversified defense cascades by LHP1
We next examined the biological processes directly influenced by LHP1-mediated H3K27me3.The scatter plots in Fig. 4a present the changes in H3K27me3 and functional gene expression in common wheat.Genes with decreased H3K27me3 levels and increased expression levels are likely directly repressed by LHP1.Consistent with the phenotype of the Talhp1-abd mutant (Fig. 2b), a group of wellcharacterized genes involved in flowering and/or floral development are targeted by LHP1, including FUL2 (from subgenomes A, B, and D), WLHS1 (from subgenome B), WSOC1 (from subgenome B), and VRN1 (from subgenome A) (Supplementary Fig. 4) 49,50 .Many LHP1 target genes encode defense-related proteins, including those involved in the perception of pathogens and signaling pathway cascades.In addition to thoroughly characterized defense genes, including the chitin sensor-encoding CERK1, which is responsive to a wide range of pathogens 51,52 , PR1, which confers resistance to a broad spectrum of pathogens 53 , and stress-responsive WRKY TF genes, the expression of the recently reported four subgenome-diversified gene clusters mediating the synthesis of defense-related metabolites 54 was apparently activated in the Talhp1-abd mutant (Supplementary Fig. 4).These results were consistent with the enriched functions among the LHP1affected genes (Fig. 4b).The enrichment of H3K27me3 in defense genes was much more apparent in wheat than in rice or Arabidopsis (Fig. 4c); most of these genes were from expanded families in Triticeae species (Fig. 1h).The expansion of defense-related gene families in Triticeae is likely associated with increased adaptability.Additionally,  the enrichment of H3K27me3 among these genes represents an epigenetic buffering system that modulates environmental plasticity.These findings imply that LHP1 preferentially represses subgenomedivergent defense pathways mediated by gene families that expanded in Triticeae.

H3K27me3 repression of subgenome-diversified defense genes is eliminated by stripe rust infection
We next focused on the type of pathogen response preferentially regulated by LHP1.We designed a statistical approach that integrates publicly available transcriptomic data.The underlying principle of this approach is that effectors in the same pathway likely trigger overlapping downstream cascades, which may be reflected by the similarity of the transcriptomic changes.For example, if mutations in two genes trigger similar differential expression patterns, it is likely that these two genes are functionally relevant.We obtained the transcriptomes of samples resistant and susceptible to various pathogens.For each pathogen, the transcriptomic changes were compared between the resistant and susceptible samples.The correlation between the transcriptomic changes revealed by these pairwise comparisons and Talhp1-abd mutant-induced changes was assessed (Fig. 5a).The sample with transcriptomic changes that were most enriched with Talhp1abd mutant-induced changes was identified as a near isogenic line (FLW29) containing a locus mediating the resistance to Puccinia striiformis f. sp.tritici (Pst) 55 , which causes stripe rust, one of the most widely destructive wheat diseases 53,56 .In other words, the lack of LHP1 triggered similar transcriptomic changes as the introduction of a Pstresistant locus.
We therefore investigated the effect of Pst infection on genomewide H3K27me3 modifications and the transcriptome.Comparing preand post-inoculation H3K27me3 levels in seedlings revealed an overall decrease in H3K27me3, which was highly consistent with H3K27me3 decrease in Talhp1-abd mutants (blue dots in Fig. 5b).Only targets with reduced H3K27me3 levels after inoculation and in Talhp1-abd mutants were primarily associated with defense responses (Fig. 5c), implying that these defense-associated H3K27me3 loci may be controlled by LHP1.By integrating Pst-induced transcriptome changes, we detected elevated transcription of a number of reported defense genes, accompanied by reduced H3K27me3 levels post-inoculation (Fig. 5d).The genomic tracks in Fig. 5e illustrates Pst infection-triggered H3K27me3 reduction and expression induction of the homologs of well-characterized defense genes.Thus, H3K27me3 represents an epigenetic buffer system that prevents stochastic processes from altering development by repressing defense genes mostly divergent across subgenomes.Pathogen infection unleashes pre-existing latent subgenomic variations by eliminating buffer activity.

Knocking out LHP1 confers resistance to wheat stripe rust
Despite the association between H3K27me3 and defense gene expression in common wheat revealed in this study and reported previously 25,26,48 , correlation does not imply causation 57 .We inoculated Talhp1-abd wheat plants with Pst race CYR32.Compared with the control, Talhp1-abd seedlings were highly resistant to Pst, with significantly fewer stripe rust sporulation sites and an apparently lower Pst DNA-to-host DNA ratio at 14 days post-infection (dpi) (Fig. 6a).The fungus was visualized by staining leaf segments with wheat germ agglutinin at 2 dpi.The average hyphal area was significantly smaller (less than half) for the Talhp1-abd plants than for the JW1 plants (Fig. 6b).Reactive oxygen species (ROS) are key signaling molecules that enable cells to rapidly respond to external stimuli 58 .We quantified the ROS content at 2 dpi via 3,3′-diaminobenzidine (DAB) staining, which showed that the H 2 O 2 level was higher in the Talhp1-abd mutant than in the control (Fig. 6c).These results indicate that the CRISPR-Cas9 editing of LHP1 leads to stripe rust resistance associated with ROS-stimulated defenses, highlighting the regulatory role of LHP1 on the response of polyploid wheat to environmental stresses.
Considered together, the results of the present study revealed the repressive effects of LHP1 on the expression of subgenome-diversified defense genes, which protect developmental processes from stochastic external changes and potentially promote neutral drift and diversity.Pathogen infections adversely affect the buffering system, enabling previously unavailable phenotypic variants to surface, thereby releasing defense cascades and facilitating the timely fixation of favorable mutations in particular environmental niches (Fig. 7).The effects of the epigenetic buffering and release of subgenome diversity lead to developmental robustness and the timely manifestation of latent phenotypes, which helps to explain the selective advantage and phenotypic plasticity of polyploid wheat.

Discussion
In this study, we clarified how the epigenetic buffering system coordinates internal genetic changes in response to external stimuli in polyploid wheat.Accordingly, we addressed the long-standing question of how genetic diversity due to polyploidization is interpreted and regulated to mediate environmental plasticity.These findings clarified the LHP1 function in terms of evolution, development, and practical applications.

LHP1 is the epigenetic anchor linking subgenome diversity and adaptive evolution in polyploids
During long-term evolution, epigenetic buffering systems may help populations reach a local fitness optimum, while also promoting neutral accumulation of potentially selectable polymorphisms in the absence of internal or external interventions.External stresses reveal cryptic phenotypes by releasing genetic variations, which are converted to a non-neutral state, ensuring that deleterious variations are periodically purged from the population and favorable variations are fixed.Because rare combinations of variations may produce a new and advantageous phenotype, increases in the diversity of allopolyploids increase the likelihood that plants will manifest advantageous traits pushing the population upward on the adaptive landscape.This represents a molecular mechanism by which adaptive peak shifts occur without adaptive valleys.Therefore, the epigenetic buffering system coupled with environmental stimuli likely shaped the trajectory of subgenome diversification and adaptive evolution in polyploid wheat.
The enrichment of H3K27me3 in subgenome-diversified genes was also observed in other polyploids, including cotton, Brassica species, and Arabidopsis, but with relatively weaker signal than those in hexaploid and tetraploid wheat species (Supplementary Fig. 5).Thus, the PcG-mediated epigenetic control of subgenome diversity is likely a common mechanism among polyploids.

LHP1 ensures appropriate developmental responses to external stimuli
The buffering effect of LHP1 on subgenomic diversity enables common wheat to resist the effects of stochastic processes and develop normally, while also contributing to defense responses to pathogen infections, ultimately resulting in environmental plasticity.These seemingly diverse effects of LHP1 are readily encompassed by a simple epigenetic framework, which presents a rich subject for further analysis.
Polyploidy is a hallmark of cancer 59 .The high plasticity of cancer cells in tumors has been attributed to high genetic heterogeneity 59 , which is similar to the close association between polyploid plasticity and subgenomic diversity.Polycomb group proteins have been implicated in cancer development and progression.Their catalytic subunit recurrently mutates in several forms of cancer and is highly expressed in numerous other types of cancer 60,61 .The oncogenic role of PcGs is under increased scrutiny because of the potential utility of these proteins for developing cancer therapeutics 60,61 .However, the underlying mechanism remains unclear and may primarily involve transcriptional repression 61 .It will be interesting to see whether our findings regarding the translation of genetic diversity to environmental plasticity and adaptive evolution via the buffering effect of PcGs are applicable to cancer cells.

Deliberate manipulation of LHP1-mediated epigenetic buffering potentially releases latent genetic diversity
The findings presented herein provide intriguing insights relevant to exploiting intrinsic immunity to enhance defenses.Identifying resistance genes or loci in wild relatives and introducing them into modern cultivars has been the main strategy used to improve common wheat defense responses.It is thought that common wheat lacks sufficient diversity to cope with environmental stresses because of polyploidy-related  bottlenecks 31 .The genomes of wheat and other crops contain a large proportion of inactive or silent genetic loci, many of which are related to key agronomic traits and stress resistance and would considerably increase genetic diversity if properly modulated 62 .In the current study, we demonstrated that although the common wheat cultivar 'JW1' is highly susceptible to stripe rust, modifying an epigenetic factor is sufficient for conferring resistance to stripe rust.The customized removal of gatekeepers provides a strategy for activating otherwise silent intrinsic defense responses.In addition, manipulating the epigenetic buffering system may facilitate the detection of key defense-related genes that are currently hidden in the wheat genome.
The introduction of an external superior allele does not necessarily result in robust superior phenotypes in all environments.Our data suggest that LHP1 is an epigenetic gatekeeper responsive to environmental changes and may ultimately determine the successful phenotypic expression of latent genetic diversity.Further characterizing the epigenetic buffering specificity may promote the adaptation of superior alleles in different hosts and various ecological niches.
Taken together, the findings of this study imply that the deliberate manipulation of the epigenetic-based storage, buffering, and occasional release of genetic variations will facilitate the detection and exploitation of useful silent loci in modern cultivars and promote the breeding of improved crops and food sustainability.

Plant materials and growth conditions
Common wheat (Triticum aestivum) variety 'JW1' seeds were surfacesterilized via a 10-min incubation in 30% H 2 O 2 and then thoroughly washed five times with distilled water.The seeds were germinated in water for 3 days at 22 °C, after which the germinated seeds with residual endosperm were transferred to soil.The seedlings were cultivated under long-day conditions.

Constructs for gene editing and wheat transformation
To modify all three TaLHP1 copies, we used two sgRNAs that target conserved regions (target1: AGGTCCTATGGCAAGCGCAA, target2: GAGCAAGCAGCAGGAGAGGT).The sgRNAs were identified by CasOT 63 .Synthesized oligos for target-specific sgRNAs were annealed and cloned into the pBUE411 vector.

Agrobacterium tumefaciens-mediated common wheat transformation
After sequencing the target sites, the binary vector was transformed into the wheat cultivar JW1 by Agrobacterium tumefaciens-mediated transformation 64,65 .Briefly, wheat spikes were collected at anthesis, harvested 14-16 days postanthesis (DPA) and sterilized with 75% ethanol for 30 s, then with 1% sodium hypochlorite (NaClO) for 15 min, and finally rinsed 5 times with sterile water under aseptic conditions.Immature embryos were isolated and incubated with Agrobacterium strain EHA 105 for 5 min.After cocultivation at 25 °C for 2 days in darkness, the embryonic axes were removed with a scalpel, and the scutella were transferred onto plates.After 5 days, the tissues were then transferred to callus selection medium for 2 weeks.The immature embryos were then placed on an induction medium for 3 weeks.The calli were then differentiated under continuous illumination (5000 lx) with fluorescent lights at 25 °C for 14 days.The regenerated shoots were transferred to the root elongation medium.The rooted plantlets were then transferred into pots and grown in growth chambers, where they were cultivated at a temperature of 20 °C with a light intensity stronger than 60,000 lx and a night temperature of 16 °C.Transgenic wheat plants were generated by Professor Genying Li (Shandong Academy of Agricultural Sciences, China).

Mutant screening
To determine the editing efficiency, we amplified the two target sites in T 0 transgenic plants for Sanger sequencing and found that the first target site was not edited, and the second target site was edited in all three subgenomes.Therefore, in the subsequent genotyping, we used Hi-TOM 66 (http://hi-tom.net/hi-tom/)to detect only the second target sites.We amplified the second target site with site-specific primers (forward: GGAGTGAGTACGGTGTGCCGTTCAGATCCTCGGTCTCT, reverse: gagttggatgctggatggATTTGAGCTGCCCTCCTGTGT, tsingke), then the products of the first-round PCR were sent to Hi-TOM sequencing.We selected two homozygous mutants (lhp1-abd-18 and lhp1-abd-19) that resulted in premature termination of the LHP1 protein for further research.The results in the main text are obtained using lhp1-abd-19.All major conclusions are consistent using the other mutant line (Supplementary Figs.6-9).All data were generated in biological duplicates (Supplementary Data 2).

Puccinia striiformis f. sp. tritici inoculations
We collected fresh urediospores of Pst race CYR32 from the leaves of wheat cultivar MingXian169.We dispersed the fresh spores in tap water and spread equal contents on the second leaf of control and transgenic wheat plants.The infected plants were incubated in a humid chamber overnight at 16 °C in the dark.We transferred the infected plants to a growth chamber (16 °C at 16 h light/8 h dark, 80% RH) to allow the growth of rust pathogens.At 48 h post-inoculation (hpi), we stained the leaf samples with 3',3'-diaiminobenzidine (DAB) (Merck Sigma-Aldrich, Shanghai, China) to examine the accumulation of H 2 O 2 .Only the site where an appressorium had formed over a stoma was considered a successful penetration.We stained the Pst hyphae with wheat germ agglutinin conjugated to Alexa-488 (Invitrogen, Waltham, MA, USA) 67 .At 14 dpi, we scanned the Pst-infected leaves to examine   7 | Model illustrating the LHP1-mediated conditional repression of subgenome-diversified defenses, which leads to developmental plasticity and facilitates adaptive evolution in allopolyploid wheat.The expression of subgenome-diversified genes is repressed by LHP1-mediated H3K27me3 under normal conditions, which protects normal development from stochastic environmental changes and promotes neutral or nearly neutral variations.Under stress conditions, the epigenetic repression is eliminated and the hidden subgenome diversification is manifested, thereby facilitating the timely fixation of favorable changes in different evolutionary niches.The extensive diversification of allopolyploid wheat increases the likelihood the population will be pushed upward on the adaptive landscape.Thus, LHP1 helps protect and potentially promote the accumulation of subgenomic diversity, while also ensuring appropriate responses to environmental stimuli, leading to developmental plasticity and facilitating adaptive evolution in polyploid wheat.

ChIP and RNA sample preparation
At the two-leaf stage, seedling leaves of 14-day Pst-inoculated and control were either frozen in liquid nitrogen for an RNA isolation step or vacuum-infiltrated with a formaldehyde cross-linking solution (0.4 M sucrose, 10 mM Tris [pH 8], 1 mM EDTA, 1% formaldehyde) for the ChIP assay.
For spike-in ChIP, chromatin from Arabidopsis was used as a spike-in control.Briefly, fragmented chromatin from all samples (including wheat and Arabidopsis) was isolated according to the normal ChIP protocol.Then, 40 µl of chromatin from all samples was reverse-crosslinked, digested with proteinase K, and subjected to DNA purification.The concentration of purified DNA was measured with Qubit dsDNA High-Sensitivity Assays (Invitrogen).Depending on the DNA concentration of each sample, chromatin from Arabidopsis and wheat are mixed together in a ratio of 1:100 for immunoprecipitation.The rest of the procedures are the same as the normal ChIP.

LASSO regression analysis
To determine which of the 191 factors would be useful for predicting subgenome-diversified genes, we used the R package glmnet 69 to perform the least absolute shrinkage and selection operator (LASSO) regression model analysis for feature selection: where Y i is the diversity for gene i, A ij is the jth factor that targets the ith gene and β j is the fitted coefficient for the jth factor.M is the total number of factors.γ is the intercept term calculated by the model.The goal of the algorithm is to minimize: where λ controls the strength of the L1 penalty.N is the total number of.LASSO analysis combined with 10-fold cross-validations was conducted to choose the most useful predictive features, and the model family we used was "binomial".

LHP1 ChIP-seq assay in wheat protoplasts
ChIP assays using JW1 leaf protoplasts were performed.JW1 plants were grown on soil under 16 h light/8 h dark conditions for 2 weeks before protoplast isolation.Approximately 30 µg of pMD19-T plasmids containing p35S:3flag-TaLHP1-B DNA were transfected into leaf protoplasts using the PEG-mediated transfection method.After incubating the protoplasts at room temperature for 48 h under dark conditions, the protoplasts were crosslinked with 1% formaldehyde in W5 solution for 10 min on ice and quenched with 32 µl 2 M glycine for 5 min.Protoplasts were collected by centrifuging at 600 × g for 2 min at 4 °C, washed with 500 µl W5 solution once and collected again.Protoplasts were lysed in 120 µl of room temperature lysis buffer (50 mM Tris-HCl pH 8.0, 10 mM EDTA, 1% (wt/vol) SDS, 1 mM PMSF, 1X protease inhibitor cocktail) by vortexing.Total lysates containing chromatin were subjected to sonication by Bioruptor until the chromatin was fragmented into 300-to 500-bp fragments.Another 400 µl RIPA ChIP buffer (10 mM Tris-HCl pH 7.5, 140 mM NaCl, 1 mM EDTA, 0.5 mM EGTA 1% (vol/vol) Triton X-100, 0.1% (wt/vol) SDS, 0.1% (wt/vol) Na deoxycholate, 1 mM PMSF, 1X protease inhibitor cocktail) was added to the lysates.The lysates were centrifuged at 12,000 × g for 10 min at 4 °C, and the supernatant was transferred to a new tube.Another 410 µl of RIPA ChIP buffer was mixed with the remaining pellet, and centrifugation, as described above, was performed again to obtain the second supernatant.The two rounds of supernatant were pooled, and the volume was brought to 1 ml with RIPA ChIP buffer.Keep 100 µl of chromatin as 10% input.Then, 20 µl agarose beads conjugated with anti-Flag antibody (A2220, Sigma) were added to the chromatin suspension and incubated for 2 h at 4 °C.After binding with chromatin, the beads were subsequently washed twice with RIPA buffer (10 mM Tris-HCl pH 7.5, 140 mM NaCl, 1 mM EDTA, 0.5 mM EGTA 1% (vol/vol) Triton X-100, 0.1% (wt/vol) SDS), LiCl buffer (0.25 M LiCl, 1% (wt/vol) Na deoxycholate, 10 mM Tris-HCl pH 8.0, 1% NP-40, 1 mM EDTA) once, and TE (10 mM Tris-HCl pH 8.0, 10 mM EDTA) buffer once.The protein-DNA complexes were eluted from beads by adding 150 µl of complete elution buffer (20 mM Tris-HCl pH 7.5, 5 mM EDTA, 50 mM NaCl, 1% (wt/vol) SDS, 50 mg/ml proteinase K) for 2 h at 68 °C with agitation at 1300 rpm.The eluate was then transferred to a new tube.The beads were eluted again with 150 µl of elution buffer (20 mM Tris-HCl pH 7.5, 5 mM EDTA, 50 mM NaCl) for 5 min.The two rounds of eluates were combined.During the elution step, the input DNA was prepared by adding 200 µl elution buffer and 7.5 µl proteinase K (20 mg/ml) and incubating at 68 °C for 2 h.ChIP DNA was extracted with phenol:chloroform (1:1), precipitated with ethanol and resuspended in TE buffer to prepare the ChIP-seq library using the

Fig. 1 |
Fig. 1 | Prediction of the core trans-factors regulating subgenomediversified genes.a Number of sc-triads (1:1:1) and non-sc-homoeologs with varying number of homologs across subgenomes in common wheat.b Workflow for predicting the core factors regulating sc-triads (1:1:1) and non-sc-homoeologs across subgenomes.c LASSO coefficients of the trans-factor binding-related regulation of non-sc-homoeologs (top) and sc-triad (bottom) genes.The coefficient value represents the relative importance of the trans-factor contributing to regulation.d Enrichment of different homoeologous groups among LHP1-targeted genes.e Fractions of homoeolog groups targeted by LHP1.f ChIP-seq read densities

Fig. 2 |
Fig. 2 | Direct repression of subgenome-diversified genes by LHP1-mediated H3K27me3.a CRISPR-editing sites in Talhp1-abd lines.The positions in gray represent the reference SNVs differentiating homoeologs.The region from 573 to 597 is the CRISPR target site.The type of mutation is listed on the right.b Representative developmental phenotypes of the Talhp1-abd mutant.Left: heading under long-day conditions, with red arrows indicating spikes; scale bar, 2 cm.Right: growth at 60 days after planting; bar, 10 cm.c H3K27me3 changes in the Talhp1-abd mutant.The x-axis represents the average read densities in H3K27me3 target genes in the wild-type and Talhp1-abd samples, whereas the yaxis represents the log 2 -transformed changes (M value) in H3K27me3 in Talhp1abd.d Enrichment of homoeologous groups among the genes with differential H3K27me3 changes (M value) in Talhp1-abd.Dot colors represent the enrichment score, with all high-confidence genes as the background.e Enrichment of genes with up-and downregulated expression among the genes with differential H3K27me3 changes [represented by M = log 2 (fold-change)], with all highconfidence genes as the background.

1 Fig. 3 |
Fig. 3 | Quantitative association between LHP1 binding and polyploidypromoted subgenome diversity.a Association between captured subgenome variations and polyploidy-promoted diversity.The captured variations are reflected by the differences between subgenomes within the same cell.The polyploidypromoted diversity was evaluated on the basis of population diversity (π).b Enrichment of the binding of trans-factors to the four regions shown in panel (a).cEnrichment of epigenetic marks in subgenome-specific regions associated with high population diversity.The epigenetic data were generated using Chinese Spring (CS) and JW1.d Genomic tracks illustrating the genomic distribution of population diversity (π), subgenome-specific regions, and the loci of LHP1-binding sites and H3K27me3.

Fig. 4 |
Fig. 4 | LHP1 preferentially represses defense cascades enhanced in wheat.a Scatter plot of the changes in the expression and H3K27me3 of functional genes in Talhp1-abd.Different colors and shapes represent different functions.b Clustering of functional domains enriched among the genes directly regulated by LHP1 with decreased H3K27me3 and increased expression in Talhp1-abd.c Enrichment of the top-ranked terms assigned to common wheat H3K27me3 target genes in rice and Arabidopsis.

2 NatureFig. 5 |
Fig. 5 | Stripe rust infection induces subgenome-diversified defense gene expression by removing H3K27me3.a Enrichment of Talhp1-abd-induced genes among the genes more highly expressed in samples resistant to various pathogens (represented by different colors) than in the corresponding susceptible samples.b Scatter plot representing H3K27me3 changes post-inoculation and in the Talhp1abd mutant.Reported defense genes are marked in red and their expression is characterized and labeled in panel (d).c Enriched functional domains among the genes with H3K27me3 levels commonly or uniquely decreased in Talhp1-abd.d Stripe rust-induced changes in the expression of genes.Reported defense genes marked in panel (b) are labeled.e Genomic tracks illustrating subgenomediversified defense gene expression and H3K27me3 changes in Talhp1-abd.

Fig.
Fig.7| Model illustrating the LHP1-mediated conditional repression of subgenome-diversified defenses, which leads to developmental plasticity and facilitates adaptive evolution in allopolyploid wheat.The expression of subgenome-diversified genes is repressed by LHP1-mediated H3K27me3 under normal conditions, which protects normal development from stochastic environmental changes and promotes neutral or nearly neutral variations.Under stress conditions, the epigenetic repression is eliminated and the hidden subgenome diversification is manifested, thereby facilitating the timely fixation of favorable changes in different evolutionary niches.The extensive diversification of allopolyploid wheat increases the likelihood the population will be pushed upward on the adaptive landscape.Thus, LHP1 helps protect and potentially promote the accumulation of subgenomic diversity, while also ensuring appropriate responses to environmental stimuli, leading to developmental plasticity and facilitating adaptive evolution in polyploid wheat.